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Abstract 

This paper presents a new method to recover the relative pose between 
two images, using three points and the vertical direction information. The 
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vertical direction can be determined in two ways: 1- using direct physical 
measurement like IMU (inertial measurement unit), 2- using vertical van- 
ishing point. This knowledge of the vertical direction solves 2 unknowns 
among the 3 parameters of the relative rotation, so that only 3 homologous 
points are requested to position a couple of images. Rewriting the copla- 
narity equations leads to a simpler solution. The remaining unknowns res- 
olution is performed by an algebraic method using Grobner bases. The ele- 
ments necessary to build a specific algebraic solver are given in this paper, 
allowing for a real-time implementation. The results on real and synthetic 
data show the efficiency of this method. 

1 Introduction 

This paper presents an efficient solution to the relative orientation problem in cali- 
bration setting. In such a situation, the intrinsic parameters of the camera, e.g. the 
focal length, the camera distortion are assumed to be a priori known. In this case 
the relative orientation linking two views is modeled by 5 unknowns: the rotation 
matrix (3 unknowns) and the translation (2 unknowns up to a scale). Its resolution 
using only five points, in a direct and fast way, has been considered as a major 
research subject since the eighties f2T\ up to now [|29l, [|20ll, [|27ll, [fT6ll. Il3l. [fl^. 
In this paper we use the knowledge of the vertical direction to solve the relative 
orientation problem for two reasons: 

1- the increased use of MEMS -IMU (inertial measurement unit) in electronic per- 
sonal devices such as smart phones, digital cameras and the low price IMU. The 
sensors fusion (camera-IMU) is not the goal of this paper, as many authors have 
shown the advantage of coupling them [[TTII . In MEMS -IMU the accuracy of head- 
ing (rotation around the vertical axis Z) is worse than for pitch (rotation around 
X axis) and roll (rotation around Y axis), due to the strength of the gravity field, 
which has no effect on a rotation around the vertical axis. Thus the new method 
presented in this paper takes a considerable benefit from a combination of data 
from MEMS-IMU and from use of 3 homologous points, that strengthen the very 
weakness of IMU data. 

2- today very performant algorithms based on image analysis are available, that 
allow to calculate the vertical direction with high accuracy. If we have only a set 
of calibrated images we can also determine the vertical direction using vanishing 
points extraction. A lot of algorithms [|2l, [[T9ll , [25|, on such topics exist in the 
literature. These algorithms are very useful in urban and man-made environments 
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[1301, m, Bn, [m. 

The use of the vertical direction so as to reduce the disparity between two frames, 
to simplify 3D vision, has already been considered by [31]. But most papers use 
a fixed stereoscopic baseline, and here we consider that we have no knowledge 
about it. Furthermore, most paper ||3T| try to solve the problem using iterative 
methods or non minimal settings (e.g. more than three points). 



2 Our contribution to the relative orientation prob- 
lem 

The main contribution of this paper is to provide an efficient algorithm to estimate 
the relative orientation using the vertical direction as an external information in 
the minimal case, using 3 points. Once the vertical direction is defined, we inject 
this information in relative orientation, based on coplanarity equation. The knowl- 
edge of the vertical direction removes 2 degrees of freedom to the problem of the 
relative orientation. Therefore it will be enough to have only 3 homologous cou- 
ples of points to solve for the 3 other unknowns: two parameters of the baseline 
because it is up to a scale and the angle of rotation around the vertical axis. These 
coplanarity constaints can be written as a system of polynomial equations. Hence, 
we solve these equations using the Grobner bases in a direct way. The possibility 
to build a solution with only 3 points is an obvious advantage in terms of compu- 
tation time, in particular when sorting the undesirable solutions by classic robust 
estimators such as Ransac (RANdom SAmple Consensus)[8]. In the Section |6] 
we show that the new 3-point method provides better accuracy and robustness to 
noise on relative orientation estimation. 

The paper is organized as follows. In the section [3] we present the geometric 
framework of our system. Section |4] rewrites the coplanarity constraint using the 
vertical direction knowledge. The resolution of polynomial system with the help 
of Grobner bases is described in Section [5| The assessment of the algorithm in 
noisy conditions is studied in Section |6.1[ where the 3-point algorithm is com- 



pared to the well known 5-point algorithm. In Section 6.2 a comparaison with 
real image database is performed. 
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Figure 1: Coordinate systems and geometry overview. The vector V^er is the 
vector of vertical vanishing point and pierces the image plane in v. Ryer is define 
in Section 1431 

3 Coordinate systems and geometry framework 

The classical coordinate system of camera (cf. figure [T]) used in computer vision 
has been chosen [11]. In this camera system {X^am^ Ycam, Zcam), the focal plane 
is Zcarn = F , F being the focal length. Given the calibration matrix K (a 3x3 
matrix that includes the information of focal length, skew of the camera, etc.), the 
view is normalized by transforming all points by the inverse of K, m = K~^m, 
in which m is a 2-coordinates point in the image. Thus the new calibration matrix 
of the view becomes the identity matrix. M is the object point. In the rest of the 
paper we suppose that all image 2D-coordinates of the point are normalized. For 
a stereo system in relative orientation, the center of the world space coordinate 
system is the optical center C of the left image, with the same directions of axes. 
The world coordinate system is denoted by (X^, Y^,, Z^). In this system the Yw 
axis is along the physical vertical of the world space. 
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4 Using the vertical direction knowledge for relative 
orientation 



4.1 Use the IMU information 

If we have of an IMU coupled with the camera, we need only to know the rotation 
angle (a) around X axis and Z axis (7) based on our coordinates system. So the 
rotation matrix equals: 



cos 7 — sin 7 
sin 7 cos 7 
1 



1 
COS a — sin a 
sin a cos a 



(1) 



4.2 Use the information given by vertical vanishing point 

If we only have a set of calibrated images of a man-made environment we can 

extract the vertical direction using vertical vanishing point. Let us suppose that 
— > 

Vyer be the vector joining C to the vanishing point in the image plane expressed in 
the camera system, and ¥^,{0 ,1 ,0) be the Y axis of the world system ((see figure 
1 1. We perform the rotation that transforms Ker into Y^. Thus, we determine the 
rotation axis 'to and the rotation angle 9 in the following way: To = V^er ® Y^, 
after simplification and normalisation 'uj = ,0 , ^^], where d = ^/V^ + 

, 9 = arccos (Ker ■ Yw), so after simplification, 9 = arccos (Vy). Using Olinde- 
Rodrigues formula we get the following rotation matrix : 

Ryer = Icos9 + sin 9[uj]^ + {1 - cos 9)uj ^u. (2) 

The rotation (Rver) given by equation [T] or |2] is then applied to all 2D points ob- 
tained in each image, m is replaced by Rver^h. 



4.3 Rewriting the coplanarity constraint 

First, we recall that for a pair of homologous points m} and iv? of a pinhole cam- 
era, the constraint on these 2 points is expressed by the equation of coplanarity: 



E 



mi 



(3) 
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where E is a 3x3 rank-2 essential matrix [fTTTl . We can also express this constraint 
by the equation |4j 
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(4) 



However, if we apply the rotation (Ryer) obtained in equation|2]to all homologous 
points, before we take in account this constraint (equation the rotation R is ex- 
pressed in a simpler way, as it remains only one parameter of rotation to estimate, 
the angle around the Y axis (vertical axis). Thus: 



Rd 



cos (f) — sin I 

1 
sin 6 cos 6 



(5) 



Using t = tan |, we replace cos by (1 — + t^) and sin by 2t/(l + 

The new coplanarity equation is rewritten as: 

i-2mlTyt + m2(T,(l - t^) + 2T,i)- 
nilTyil-t'))ml + {ml{l + t')T,+ 



ml{l + t')TMl + {mlTy{lt') + 
mmr^t - T,{1 - t^)) - 2mlTyt)ml = 0. 



(6) 



3 pairs of homologous points allows for instancing equation 6 as {/2, /s, U} 
with remaining unknowns Tx,Ty, Tz and t. The corresponding base is only com- 
posed from two degree of freedom since no scale modeling has been yet per- 
formed. Therefore it is necessary either to fix a component of the base to 1, either 
to add the constraint of normality. We choose this last one: /i = T| + + T"^ — 
1 = 0. The advantage is that it allows to get a more general modeling. We have 
therefore a system of 4 polynomial equations of degree 3 {fi, f2, fs, fi}- Now we 
describe the direct resolution of this polynomial system using the Grdbner bases. 



5 Resolution of the relative orientation equation us- 
ing Grobner bases 

We recall first the basic definitions of Grobner bases, and also the link between 
Grobner bases and linear algebra. Then, we use these concepts to derive a specific 
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algorithm to compute the Grobner basis of the system of polynomials defined in 
Section O 

5.1 Properties of Grobner basis 

The notion of Grobner basis was introduced by B. Buchberger, who gave the first 
algorithm to compute it (see [4J). This algorithm is implemented in most gen- 
eral computer algebra systems like Maple, Mathematica, Singular [jTOl . 
MACAULAY2 [9], Cocoa [5J and Salsa software [22j. Let R = K[xi, . . .,Xn] 
be a polynomial ring where K is an arbitrary field. Let fi, . . . , fk G i? be a 
sequence of k polynomials and let / = (/i, . . . , /fc) be an ideal of R gener- 
ated by the //s. We need also a monomial ordering on R. We recall here 
the definition of the degree reverse lexicographic ordering (DRL), denoted by 
-<, which is an especial monomial ordering having some interesting computa- 
tional properties. For this we denote respectively by deg(m) (resp. degj(m)) 
the total degree (resp. the degree in x,) of a monomial m. If m and m' are 
monomials, then m ^ m' if and only if the last non zero entry in the sequence 
(deg^(m') — deg]^(m), . . . , deg„(m') — deg„(m), deg(m) — deg(m')) is negative 
(see L7J). 

Let in(/) G i? be the initial (greatest) monomial of a polynomial f E R with 
respect to -< and in(/) = (in(/) | / G /) be the initial ideal of /. 

Definition 5.1 (Grobner basis) A finite subset G C I is a Grobner basis of I 
w.r.t. -< if{m{G)) = in(/). 

Definition 5.2 (Reduced Grobner basis) A Grobner basis G of I is called re- 
duced if for all g E G, g is monic and no monomial ofg lies in {m{G \ {g})). 

Proposition 5.1 ([TJ, Proposition 6, page 92) Every ideal has a unique reduced 
Grobner basis. 

5.2 Macaulay matrix 

We recall now the definition of a Macaulay matrix and we explain who we could 
use it to compute the Grobner basis of an ideal. With the notations of above 
subsection, we consider the ideal / generated by the //s and ^ be DRL monomial 
ordering. We suppose that we know the maximum degree d of monomials which 
appear in the representation of the elements of the Grobner basis of / in terms of 
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the /j's (in Subsection |5.3[ we show how to compute such a degree for the ideal 



generated by polynomials defined in Subsection 4.3 1. Note that this degree is the 
maximum degree of monomials which appear in the computation of the Grobner 
basis of /. 

We can build the Macaulay matrix Mrf(/i, . . . , f^) (for short we denote it by 
Mrf) as follows: Write down horizontally all the monomials of degree at most d, 
ordered following -< (the first one being the largest one). Hence, each column of 
the matrix is indexed by a monomial of degree at most d. Multiply each fi from 1 
to k by any monomial m of degree at most d — deg(/j), and write the coefficients 
of mfi under their corresponding monomials, thus giving a row of the matrix. The 
rows are ordered: row mfi is before ufj if either i < j or i = j and m ^ u. 

monomials of degree at most d 
: / \ 

\ I 

For any row in the matrix, consider the monomial indexing the first non-zero 
column of this row. It is called the leading monomial of the row, and is the leading 
monomial of the corresponding polynomial. 

Gaussian elimination applied on this matrix leads to a Grobner basis of I (see 
[fTSlP . Indeed, call the Gaussian elimination form of M^, such that the only 
elementary operation allowed for one row is the addition of a linear combination 
of the previous rows. Now, consider all the polynomials corresponding to a row 
whose leading term is not the same in Md and Md, then the set of these polyno- 
mials is a Grobner basis of /. 



5.3 Constructing the specific Macaulay matrix 

In this subsection we describe a general algorithm to compute the Grobner basis of 



the system of polynomials defined in Subsection 4.3 It is worth noting that when 
the coordinates of the input points change, only the coefficients of polynomials 
change. Thus, using Lazard's approach (see the above subsection), we build a 
Macaulay matrix (and we may compute it directly when the coordinates of the 
input points change), and a Gaussian elimination on this matrix gives the Grobner 
basis of the ideal. 

Let fi, ■ ■ ■ , fi E C[Tx,Ty,Tz,t] he the system of polynomials as defined in 



Subsection 4.3 Let / = {fi, ■ ■ ■ , fi)- Our first challenge is to choose a good 
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monomial ordering. From a good monomial ordering, we mean an ordering for 
which the maximum reached degree in Grobner basis computation is minimum. 
Or in terms of complexity, we look for an ordering for which the computation has 
the optimal complexity. We choose DRL ordering because it typically provides for 
the fastest Grobner basis computations. Let us consider DRL(Tj., T^, T^, t). We 
compute first the maximum degree of monomials which appear in the computation 
of the Grobner basis of / w.r.t. this ordering. We use this degree to study the 
complexity of computing Grobner basis and also to construct the Macaulay matrix 
of / to compute its Grobner basis. For this, we homogenize the //s w.r.t. an 
auxiliary variable h and we compute the Grobner basis of the homogenized system 
for DRL(Ta;, Ty, Tz,t, h). The maximum degree of the elements of this basis is 6 
and therefore the maximum degree of monomials which appear in the computation 
of the Grobner basis of / will be 6 (see [T5] for more details). We have tested some 
other monomial orderings, and it seems that this ordering is the best one. 

Our second challenge is to build M6(/i, . . . , f^), say M. To compute such a 
matrix, we have to find the products mfi, such that a Gaussian elimination on the 
matrix representation of these products leads us to the Grobner basis of I. For 
this, we use the maximum reached degree in Grobner basis computation which 
is 6. We consider all products mfi where m is a monomial of degree at most 
6 — deg(/i). This gives 175 polynomials. Among them, there are some products 
which are useful to build M. Using the following programme in Maple, we 
could choose the useful ones: 

L:=NULL: 
AA : =A : 

for i from 1 to nops (A) do 
unassign ( ' p' ) : 
X : =AA : 

member (A [i], AA, 'p'): 
AA:=subsop (p=NULL, AA) : 
if IsGrobner (Macaulay (AA) ) then 
L : =L, i : 

else 

AA : =X : 

f i : 

od : 

where IsGrobner is a programme to test whether a set of polynomials is a 
Grobner basis for / or not, and Macaulay is a programme which performs a 
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Gaussian elimination on the matrix representation of a set of polynomials. This 
gives 65 polynomials of degree at most 6. In this case, M has a size 65 x 77. Here 
is the list of 65 polynomials which were found by this way. 

Ja, tfi, Tzfi, Tyf4, TxJa, tTzf4, tTyfi, tTxfi, 

Ti rji r rji rji _r rjil n rji rji n rTi2 £ irji rji r 
zJ-yjA-, J^zJ^xjA-, J-yJAi J-yJ-xjA-, 1 ^ J A-,I^J- zJ- yj A, 

tTzTxfA, tTyfi, tTyTxfA, tT^fA, /s, tfs, T^fs, 
TyfsjTxfs, tTzfs, tTyfs, tTxfs, T^Tyf^, T^T^fs, 
'^yfs: TyT^fs, T^fs, tT^T^fs, tTyf^, tTyT^fs, 

fs, f2, tf2, Tzf2, Tyf2, Txf2, tTzf2, tTyf2, 

tTxf2,TzTyf2, TzTxf2, Tyf2, TyTxf2, 

tTzTxf2, tTyf2, tTyTxf2, tT^f2, fi, tfi, Tzfi, 

Tyfi, Txfi, t^/i, tTyfi, tTxfi,t^fi,t^Tyfi, 

t^T^fii t^Tyfi, t^T^fi 

Remark that IsGrobner and Macaulay were written in Maple and the 
former does not use Buchberger's criterion to test whether or not a set of polyno- 
mials is a Grobner basis or not, because using this criterion is very time-consuming. 
In fact, we have used the properties that we can compute in(/) and a set of polyno- 
mials G C / is a Grobner basis for / if in(G) = in(/). This makes IsGrobner 
very fast and efficient, and allows to do the above choice in real time. 



5.4 Constructing the specific algebraic solver 

In this subsection,, we recall briefly an algebraic solver which uses a Grobner 



basis to find the solutions of the system defined in Subsection 4.3 



Thanks to the property that the division by the ideal / is well defined when we 
do it w.r.t a Grobner basis of /, we can consider the space of all remainders on 
division by / (see [7]). This space is called the quotient ring of /, and we denote 
it by A = C[Tx, Ty, Tz, t]/I. It is well-known that if / is radical then the system 
/i = ■ ■ • = /4 = has a finite number of solutions if the dimension of A as an 
C-vector space is (see LZJ, Proposition 8 page 235). We can easily check by the 



10 



function I sRadical of MAPLE that / is radical. A basis for A as a vector space 
is obtained from in(/) by (0, Theorem 6, page 234) 

B = {m I m is a monomial and m ^ in(/)} 

From computing a Grobner basis of /, we could compute in(/), which is equal to 
in(/) = (T^, Ty, T2, f) and thus the set 

B •[l,t,t ,t ,t ,t ^TziT^t^T^i ^T^t ^Tzi yT^t } 

is a basis for A as an C-vector space. Therefore, we can conclude that the system 
/i = ■ ■ ■ = /4 = has 12 solutions. Note that we have obtained these results 
for an especial coordinates of input points. We can discuss mathematically the 
correctness of these results for any set of points. But, that is out of the subject of 
this paper and the scope of this conference. We recall here briefly the eigenvalue 
method that we have used to solve the system /i = ■ • ■ = /4 = 0, see [6J, page 56 
for more details. For any / G C[Tx, Ty, Tz, t] let us denote by [/] the coset of / 
in A. We define irif : A — > A by the following rule: 

mfi[g]) = [f].[g] = {f9]eA 

Since, the ideal generated by the /j's is zero-dimensional, then A is a finite di- 
mensional C-vector space, and we can present m/ by a matrix which is called the 
action matrix of /. For any i, if we set / = Xj, then the eigenvalues of ma;, are the 
Xj-coordinates of the solutions of the system. Using these eigenvalues for each i, 
and a test to verify whether or not a selection n-tuple of these eigenvalues vanishes 
the /i's, we could find the solutions of the system. A more efficient way is to use 
eigenvectors. Let / be a generic linear form in A, then we could read directly all 
solutions of the system from the right eigenvectors of mj, see [l6l, page 64. 



5.5 Computation of final relative orientation 

After the resolution of the polynomial system, and the obtention of the parameters 
Tx ,Ty , Tz and t, it is possible to compute the finale relative orientation between 
the images. If we suppose that Rveri is the rotation matrix defined in the section 



4.2 for the image 1, and Ryer2 same for the image 2, and R^ the rotation matrix 
defined by t (equation |5]), the final relative orientation between the images 1 and 
2 is: 

R final Rver2 RdiRverly 

(7) 



Tfinai = Rver2^ T, whcre T = [Tx,Ty,Tz\K 
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6 Experiments 



The accuracy of the relative orientation resolution, using a vertical vanishing point 
and 3 tie points, is based on three factors : 

1- the accuracy of the polynomial resolution of the translation parameters (Tx , Ty , Tz), 
and of the rotation around the Y axis using the Grobner bases, 

2- the geometric accuracy for the estimation of the vertical direction, 

3- the accuracy of the algorithm on tie points in presence of noise. 

In order to evaluate the different impacts, we have in a first time worked on syn- 



thetic data in Section 6.1 then we have used real data in Section 6.2 



6.1 Performance Under Noise 

In this section, the performance of the 3 points method in noisy conditions has 
been studied and compared to the 5 points algorithm fTT\ using the software pro- 
vided by authors [|26ll . The employed experimental setup is similar to ||20ll . The 
distance to the scene volume is used as the unit of measure, the baseline length 
being 0.3. The standard deviation of the noise is expressed in pixels of a 352x288 
image as cr = 1.0. The field of view is equal to 45 degrees. The depth varies be- 
tween to 2. Two different translation values have been treated, one in X (sideway 
motion) and one in Z (forward motion). The experiments involve 2500 random 
samples trials of point correspondences. For each trial, we determinate the an- 
gle between estimated baseline and true baseline vector. This angle is called here 
translational error, and expressed in degrees. For the error estimation on the rota- 
tion matrix, the angle of {Rjrue^esUmate) is calculated, and the mean value for the 
2500 random trials for each noise level is displayed. From Figure [2[|3}|4] and [5| we 
see that the 3-point algorithm is more robust to error caused by noise in sideway 
and forward motion for estimation of rotation and translation. 

Now let us compare 3-point and five-point algorithm on a planar scene. In 
this configuration all the points of the scene in the world have the same Z (here 
equal to 2). The results for the estimation of the rotation (Figure |6]) show that the 
two algorithms provide a good determination of the rotation, but the 3-point gives 
much better results than the 5-point one for the base determination in sideway 
motion (Figure |7]). This weakness of the 5-point algorithm in planar scene has 
been discussed in [|24ll . 
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Figure 3: Error on the baseline orientation (in degrees, sideway motion). 



6.1.1 Impact of the accuracy of the vertical direction on the estimation of 
relative orientation 

We have introduced an error of to 0.5 ° on the angular accuracy of the vertical 
direction. Today for example, a low-cost inertial sensor such as Xsens-MTi [fT2ll 
gives a precison around 0.5 ° on rotation angle around X axis and Z axis (the ver- 
tical direction being Y axis). Of course, some high accuracy IMU are available, 
they may reach an accuracy better than 0.01 ° on the orientation angles if properly 
coupled with other sensors (e.g. GPS). Using an automatic vanishing point detec- 
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Figure 4: Error on the rotation (in degrees, forward motion) 
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Figure 5: Error on the baseline orientation (in degrees, forward motion). 



tion specially in urban scene, we get a very precise vertical direction (better than 
0.001 °), as it will be shown later. We have checked the impact of this accuracy on 
the determination of the rotation and the base. (Figure [10] and Figure [TT]). 



6.2 Real Example 

So as to provide a numerical example on real images, we have chosen to work on 
the 9-images sequence "entry-PlO" of the online database [|28l . In this database 
we know all the intrinsec and external parameters. First, we extracted the vanish- 
ing points on each image. We used the algorithm of [T3| because beyond its high 
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Figure 6: Error on the rotation (in degrees) in planar configuration (sideway mo- 
tion) 
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Figure 7: Error on the base orientation (in degrees) in planar configuration (side- 
way motion) 



speed, it allows an error propagation on the vanishing points according to the error 
on the segments detection. We express this error in an angular manner. The results 
of the angular errors are shown in the table 1 . As one can see it, the determination 
of the vertical vanishing point is very precise and according to the Figure [10] and 



11 it induced an error close to zero. Then, we have computed the relative orien- 



tation for 3 successive images (each time, 2 following couples of images). The 
interest points are extracted using SIFT [18] algorithm. The results are presented 
in the Figure [12} The mean value of angular errors on the rotation amounts to 0.82 
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Figure 9: Error on the base orientation (in degrees) in planar configuration (for- 
ward motion) 




(a) (b) 

Figure 10: Impact of the geometric accuracy of the vertical direction on the esti- 
mation of a) the rotation (in degrees), and b) the base orientation (in degrees) in 
sideway motion. 
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(b) 



Figure 1 1 : Impact of the geometric accuracy of the vertical direction on the esti- 
mation of a) the rotation (in degrees), and b) the base orientation (in degrees) in 
forward motion. 



Image 


Angular error on vertical direction in degree 


0000 


0.002569 


0001 


0.0066 


0002 


0.001584 


0003 


0.001443 


0004 


0.000899 


0005 


0.00115 


0006 


0.001445 


0007 


0.005018 


0008 


0.002424 


0009 


0.002223 



Table 1: Results. Vertical direction detection using the vertical vanishing point. 

degree. For the estimation of the translation, this error amounts to 1.33 degree. 
These results show clearly the efficiency and robustness of the method. 

6.3 Time Perfomance 

The resolution of the polynomial system and detection of vanishing point was 
written in C ++. With a 1.60 GHz PC the time of each resolution is about 2 fis, al- 
lowing real-time application. We may note that the selection process using RanSac 
dHl among the SIFT points is running considerably faster on 3-point than on 5- 
point algorithm. 
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Figure 12: Result on "entry-PlO" sequence. Each cell contiens the error on rota- 
tion in degrees (upper left) and error on the translation in degrees (bottom right). 

7 Summary and Conclusions 

Today, more and more low-cost personal devices include MEMS-IMU in comple- 
ment to cameras, these devices allow to provide very easily the direction of the 
vertical in the image. Furthermore, image based automatic extraction of the ver- 
tical vanishing point offers a very high accuracy alternative, if needed. So, here, 
we have demonstrated the advantage of using the vertical direction, and an effi- 
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cient algorithm for solving the relative orientation problem with this information 
has been presented. In addition to a considerable acceleration, compared with the 
classical 5 point solution, our algorithm provide a noticeable accuracy improve- 
ment for the baseline estimation. Another interesting feature improvement has 
been demonstrated: the planar scenes raise no more problem in baseline estima- 
tion. This advantageous result is due to an appropriate problem formulation using 
in a explicit way the significant parameters of the relative orientation (parameters 
of the rotation and the translation). 
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